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Abstract 

Vacuum expectation values of products of local bilinears ipip are com- 
puted in massless QED 2 at finite density. It is shown that chiral con- 
densates exhibit an oscillatory inhomogeneous behaviour depending on 
the chemical potential. The use of a path-integral approach clarifies 
the connection of this phenomenon with the topological structure of the 
theory. 
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I. INTRODUCTION 



Fermion condensates play an important role in particle physics and cosmology. In par- 
ticular, they are relevant in connection with chiral symmetry breaking 0, a phenomenon 
related to the structure of the QCD vacuum . 

Recently, it has been found in the large Nq limit Q, that the order parameter for chiral 
symmetry, < ipip >, is at high fermion densities inhomogeneous and anisotropic so that the 
ground state of quark matter has the structure of a standing wave with respect to the order 
parameter. 

At present, it is well understood that the fermion condensate < tjj(x)i/j(x) > picks contri- 
butions from non-trivial topological sectors. Then, in order to decide whether the breakdown 
of chiral symmetry occurs, one should consider instanton effects on the fermion condensate. 

Two-dimensional models like QED 2 (the Schwinger model) and QCD 2 provide a natural 
ground to study these phenomena since, although simplified, the basic features (chiral sym- 
metry breaking, non-trivial topological sectors, etc) are still present and exact calculations 
can be in many cases performed. In this context, we calculate in the present work vacuum 
expectation values of products of local bilinears, ifii/j, at finite density for the Schwinger 
model. We employ a path-integral approach which leads in a very simple way to exact 
results and has shown to be very adequate for studying non-Abelian extensions [Q- [|7| . Our 
conclusions go in the same direction as those in || and extend the results obtained in ||. 
In particular we are able to show that multipoint chiral condensates exhibit an oscillatory 
inhomogenous behavior depending on the chemical potential, as it is the case in the 1/N C 
results for QCD^ ||. The lack of translation invariance manifests through a dependence 
on differences as well as sums of spacial coordinates. This last result corrects the oscilla- 
tory behaviour computed approximately in || for the two-point correlator in the Schwinger 
model. 

Let us start by observing that Quantum Field Theories at finite fermion density can be 
studied by introducing a chemical potential leading to a quantum theory in the presence of 
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a classical background charge distribution 0- [|Tl|]. Concerning topological contributions, 
Bardakci and Crescimanno [0 proposed a natural way to take into account topologically 
non-trivial configurations, which is suitable for the study of the fermion condensates within 
the path-integral formulation. In this approach, one decomposes a given gauge field belong- 
ing to the n th topological sector in the form 

A,(x)=A(;\x)+a,(x) (1) 

with A^> a fixed configuration carrying the (whole) topological charge n and the path- 
integral variable, which accounts for quantum fluctuations and belongs to the trivial topo- 
logical sector. In this way, calculations easily workable in the n = background (as the 
evaluation of Fujikawa jacobians after a fermionic change of variables using techniques re- 
quiring a compact manifold) can be handled without problems and, at the same time, the 
contribution of topologically non-trivial sectors is properly taken into account. 



II. THE MODEL 

We start from the two dimensional (Euclidean) Lagrangian for QED 2 

C = j)* -±F%,, (2) 

where our 7 matrices are taken as the Pauli matrices, 7o = 7i = —02 so that 7^75 = 

The partition function Z for the model is 

Z = J DipD^DA^expi- J d 2 xC) (3) 

In order to include fermion density effects, it will be convenient to add still another vector 
field A b ^ which will allow for the introduction of the chemical potential in a simple way. This 
field describes an external charge density acting on the quantum system. Whenever A h ^ is 
taken as i times the chemical potential [0, it will represent a uniform charge background. 
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If one first considers a finite length (2L) distribution and then takes the L — > oo limit, 
translation symmetry breaking becomes apparent and at the same time ambiguities in the 



definition of the finite density theory are avoided (see ref. flTOl ). 



Now we can write the gauge field in Lagrangian (H) in the form 

A,(x) = A^(x) + A b l (x)+a,(x) (4) 
Hence, the Lagrangian in the n th topological sector reads 

(^) 2 + + 2i#>i£,) + £ c (5) 

In eq.(|5|) we have added a counterterm C c which will be conveniently chosen to cancel out a 
divergency arising when the background A b ^ is taken as a constant to introduce the chemical 
potential (see below). 

The partition function now reads 

Z = J2 I DipD^Da^expi- f d 2 xC {n) ) (6) 

where we have written the path-integral as a sum over all topological sectors. 

As it is well-known, in two dimensions fermions can be completely decoupled from gauge 
fields by an appropriate chiral and gauge transformation of the fermion fields. At the 
quantum level, the corresponding change in the fermionic variables is accompanied by a 
non-trivial Fujikawa jacobian fl3f . As it can be inferred from the connection between the 



chiral anomaly (and the index theorem for the Dirac operator) with the change in the 
fermionic path-integral measure, calculation of the jacobian can be more easily performed in 
compact space-time manifolds. The latter implies that only decouplings from topologically 
trivial gauge fields should be considered. With this in mind we perform the following change 
in the fermionic variables 

if) = exp(7 5 + «7 + 7 5 6 + ir] b ) \ 

4> = X exp(7 5 -irj + 7 5 b - ir] b ) (7) 



according to a standard decomposition of the vector fields 

<V = - e ^du4> + dp?) (8) 

and 

4 = -e^dvfi + d,r, b (9) 

With this change of variables the only field that remains coupled to the fermions is Aft\ 
a classical configuration carrying all the topological charge. In terms of the new variables, 
Lagrangian ([|) becomes 

£<»> = m + 4 (n) )x _ -L( F (;)) 2 - ( (D0 )2 + e^Ucf) + £ c (10) 

where we have already taken A b ^ as a constant, in order to introduce the chemical potential. 

Eq. ( |T0| ) describes the effect of the change of variables merely at the classical level. At the 
quantum level, within the path-integral approach, one has to take into account the change 
in the fermionic measure, i.e., the Fujikawa jacobian associated with transformations (^). 
The evaluation of this jacobian is standard and we just quote the result ||] 

J = exp (— J d 2 x 0D(0 + 20 (n) )) x exp (J- J d 2 x 6 D(0 + 20 (n) )) (11) 

where (p^ is accordingly defined by 

4") = -e^d^ (12) 

Concerning the Jacobian associated with the change of the bosonic variables (eqs.([|),(§)), 

Da^ = J bos D(f)Dr] , (13) 

it can be identified with the Faddeev-Popov determinant for the rj = (Lorentz) gauge, 
Jbos = detn, and it will be ignored in what follows. 

To relate with the chemical potential /z, note that fixing the fermion number density 
through a term 
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corresponds to the choice 



or, equivalent ly, 



£chem = -Z^V'ToV' 



A\ = -itfrt (15) 



ifixi rjb = (16) 



Notice that, as has been thoroughly analysed by Actor |TI[] , does not correspond to a 
pure gauge. Were it not so, the introduction of a chemical potential would not have physical 
consequences. For the same reason, one cannot gauge away C c hem from the lagrangian by 
means of the alternative choice: rji, = —ifix^, <pb = 0. In fact, this transformation would 
correspond to an unbounded (in the temporal axis) gauge one. Although transformation ([16]) 
is unbounded as well, as mentioned above one can handle this problem as it is usually done: 
puting the system in a spatial box, introducing counterterms and taking the infinite limit for 
its length at the end of the calculations. In this way the usual divergency associated to the 
electromagnetic energy carried by fermions can be eliminated by an appropiate counterterm 
|T0| . In our approach the divergency manifests through the term (p^Ucf)^ in eq . (|Tl~D . As 
stated above, this infinite contribution is cancelled out by an appropriate choice of C c . 
This counterterm is the Lagrangian counterpart of the usually employed in the Hamiltonian 
approach to handle this problem ||10]. In the canonical QFT this is tantamount to a 
redefinition of creation and annihilation operators which is equivalent to a shift in the scale 
used to measure excitations. 

Putting all this together, the partition function of the model, can be written in the form 

Z = ATE / DxD X D<t> exp(-£$) (17) 

where is the effective action in each topological sector, 

Sif f = Jd 2 x ^ + /")) X -^|rf 2 x((D0) 2 + e^)n0) + 

~h I d2 < F ^ 2 -k^\ Sx 0D(0 + 20(n)) (18) 



Note that with the choice of C c discussed above, the effective action written in terms 
of the "decoupled" fermions does not depend on the chemical potential /i. Nevertheless, \x 
reappears when computing correlation functions of fermion fields, once ip and ip are written 
in terms of the decoupled fields \ an d X through eq.(^). 

Of course, the fermionic integral in eq. flT7|) is the determinant of the Dirac operator in 



the background of a gauge field carrying topological charge n. Now, as it is well-known [12 



the Dirac operator has, for n > (n < 0), n positive (negative) chirality zero modes so that 
actually none but the n = sector does contribute to Z. However, in computing v.e.v.'s of 
products of fermion bilinears, Grassman coefficients accompanying zero-modes render non- 
trivial certain path-integrals in a given topological sector. This happens according to the 
number of bilinears appearing in the v.e.v., and then only this sector will contribute (see 
next section). 



III. CORRELATION FUNCTIONS OF FERMION BILINEARS 

We are now ready to study the behavior of chiral condensates and their dependence on 
the chemical potential /i. To this end, let us define the chiral charge changing correlators 

s + (w) = $ + if>+(w) (19) 

s_{w) = $-ip-{w) (20) 

where ip + and tp_ are the right-handed and left-handed components of the Dirac spinors. 
Thus, the fermion condensate < ipif}(w) > is the sum of the v.e.v.'s of the composites defined 
in eqs.p)-©. 

Notice that this correlator could be naively expected to vanish as a direct consequence 
of the chiral invariance of massless QED 2 . However, as explained above, when one is to 
compute < s + (w) > (< S-(w) >) using the partition function given by eqs.(p~7])-([T8|), a non 
zero value is found. In fact, only the n = 1 (n = — 1) sector will contribute to the sum over 
topologically non trivial sectors. One then has 
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< s±(tw) > = — / DxDxDcj) X±{w)x±{w) exp (±2(0(iw) + 0bO))) 

(21) 

where Z , the only non-vanishing contribution to Z, arises from the n = sector. After 
some algebra, eq.(^) becomes 

< s±(w) > = -J- J DxDx X±{w)x±{w) exp(- j d 2 x x(ift + 4 (±1) )x) 
±- J D<f ) exp(±2(<f ) (w)+t f iw 1 )) x exp (-_L(F / (± 1 )) 2 ) 
exp(^- J d 2 x(<pn<p + 2<j>e^A^) (22) 

with 

Z F = detifi (23) 

and 

Z B = det- 1 / 2 (□(□ - |^ (24) 

It can be easily proved in the /i = case that < s + >+7° an d < s - >-7° coincide. This 
in turn implies that 

< > At= °= 2 < s+ >^7° (25) 

Now, from eq.(^) we can see that 

< s± >±i= exp(±2i/xwi) < s± >±i° , (26) 

hence, one finally gets 

< ^(w) >^°= cos(2/iWi) < > M=0 (27) 



In order to have an explicit formula for eq. (27) , let us recall that 



where 7 is the Euler constant and m = e 2 /ir is the mass of the effective boson jnj]. Eq. 
coincide with that presented in || using operator bosonization rules. Appart from the 
simplicity of our derivation in the path-integral framework, it should be stressed that it is 
within this approach that the role of different topological sectors becomes apparent. This 
can be put in evidence in the calculation of N point correlators which will be discussed in 
the remaining of this section. 

According to the previous discussion, we consider the contribution of the n th topological 
sector to the N point correlation function 

N 1 N 

< IKK) >fn = -w- \ D X D X D<Pl[xP±x(w t ) x 
i=i ^0 J i= i 

N 

exp(±2X;e(«; 4 ) - Sffi) (29) 
1=1 

Here < > ±n means that the v.e.v. is computed in the in^-topological sector; P + (P_) is 
the projector on the right-handed (left-handed) subspace 

P±=\{l±ls) (30) 

and we have defined 

6{w) = <j)(w) + (j) b (w) (31) 

Now, due to Grassman integration rules, it can be easily proved that the only non vanishing 
contribution to these correlators may arise when the number of insertions equals the absolute 
value of the topological charge, ie: \n\ — N [IS]. Performing the decoupling of fermions as 



before we finally get the following expression 



N 



Here 



< II s ±( w i) >±n = F ±N (w 1 , w 2 , ...,w N )B± N (w 1 , w 2 , ...,w N ) (32) 



. N 

F ±N = Zp l J D X D X II X±X±(m) exp(- / d 2 x X (z0 + 4- {±N) )x) (33) 

i=l 



and 



B ±N = Z B l J D(f) exp(±2£0(O)exp(— J d 2 x<pn(<p + 20< ±JV >)) x 

exp (J^ J d 2 x{{U<pf + 2D0( ±7V )D0 + (D0 (±iV) ) 2 ) (34) 



We thus see that the dependence of the correlator on the chemical potential factorizes, 
resulting in 



N 



<n 



S±{Wi) > ±N 



i=l 



N N 
i=l i=l 



Wi) >% N 



(35) 



where 



N 



<Y[ s± i 



Wi) >± N = 



-me 1 /47r) exp(— 2 V] i^ (^|wi — 



(36) 



are the general minimal correlation functions for /i = [15 , 16 



We are now ready to go further and study the so-called non-minimal |]12[ , |16| , p77 



corre- 



lation functions which will be needed in order to compute multipoint vacuum condensates. 
We start considering the complete two point composite 



< ipip(x)il)i/;(y) > = < s + (x)s + (y) > + < s„(x)s + (y) > 
+ < s + (x)s_(y) > + < s_(x)s_(y) > 



(37) 



In the n = topological sector the unique contribution to eq.(^) comes precisely from the 
simplest non-minimal v.e.v. which after fermion decoupling becomes 

< s„(x)s+{y) + < s+(x)s_(y) >^°= 
(exp (— 2£/i(ari -j/i))x < X+X+(^)X-X-(y) >f + 

exp(2^! < >f ) e/^^+i^*)) (38) 



This gives a contribution to < ipip(x)ipip(y) > in the n = topological sector, of the 
form 



< ijjij(x)^(y) >^o= cos(2/x(xi - yi)) < xpip(x)ipip(y) >£J 



At=0 



(39) 
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There are two remaining contributions arising from the n = ±2 sectors. They follow directly 
from eq. (|35|) 

< s + (x)s+(y) >£f° = exp(-2i/x(xi + yi)) < s + (x)s + (y) >£f 

< s-{x)s-(y) = exp(2z / u(x 1 + yi )) < (40) 

Again 

< s + (x)s+(y) >l= 2 °=< s4x)s4y) >^°= ^ < ^(x)^(y) >g° 2 (41) 
so that we finally have 

< ^i){x)i)il){y) >^°= cos(2/i(xi - yi))x < ^(x)tl;ip(y) >%Zo + 
cos(2 / u(x 1 + Vl )) x < $i/>(x)W(y) >fj° 2 (42) 

At this point some remarks are in order: 

(i) Once again the topological structure of the theory manifests, making possible the discrim- 
ination of the contributions to fermion bilinear correlation functions from each topological 
sector. 

(ii) Eq.(|42|) exhausts all topological contributions. 

(iii) The lack of translation invariance (which is broken by the background charge distribu- 



tion) becomes apparent particularly through the last term in the r.h.s. of eq. (f42|) , which 
depends on the combination x\ + y\ 

(iv) No clustering ansatz has been needed in order to obtain these results. 

Let us finally note that if we put y\ = in eq.(|42|) we get a compact formula analogous 
to eq.flSTD, 

< i)ij){x)i)^{y) >^°= cos(2//xi) < i)^{x)i)ijj{y) > M=0 (43) 

This last result can be seen to coincide with that obtained in ref. || using cluster decom- 
position. In fact, the result reported in || corresponds just to the trivial topological sector 
and does not reproduce the contribution of n ^ sectors whenever x\ ^y\. 
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Correlators of a larger number of bilinears can be very simply obtained following the 
same procedure as above. As an example one gets for the 3-point correlator, 

n 

2 cos(2/j l (x 1 + Vi + zi)) < s + (x)s + (y)s + (z) >£l° + 
2cos(2/i(xi +y 1 -z x )) < s + {x)s + {y)s_(z) >£zj + 
2cos(2/i(xi -7/1 + 2:1)) < s + (x)s-(y)s + (z) + 

2cos(2 / u(-xi +yi + z x )) < s_(x)s+(y)s + (z) >£l? (44) 



(We have emphasized that the l.h.s. of eq.(f4"4"|) exhausts all the topological contributions by 
explicitely showing the sum over n). 

From the examples above one can easily infer the structure for the general N point 
correlator 

< ^{w v )^{w 2 ) . ..#(/) >^°= 
£cos(2ME^i)) < s+M'+M • • • s + {w N ) >TJ N + 



k ij^k 



(w k - l )s(w k ) 



s + (w k+1 ) . . . s + (w N ) +Ecos(2/,( £ w\ -w k - w{)) < s + {w l ) 

s + (w 2 ) . . . s + (w k - 1 ) 8 -(w k )8 + (w k+1 ) . . . s + (^'- 1 )s_(^') 

s + (w> +1 )...s + (w N )>»Z° N _ 2 +... (45) 

A compact expression for the non-minimal correlation functions appearing in the last 
equation is 

r s 

< II s +( w i) II S -K) >r-s°= (-me747r)"e X p(-2$> i e i Ar Mw i - *%|)) (46) 

8=1 j = l i>j 

where r + s = N (for the details see |]16| ). 

In summary, we have presented the correlation functions of fermion bilinears in QED 2 
at finite density, using a path-integral approach particularly adequate for identifying con- 
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tributions arising from different topological sectors. We have been able to exactly compute 
correlation functions for an arbitrary number of bilinears, showing its dependence with the 
chemical potential. One of our motivations was a recent work by Deryagin, Grigoriev and 
Rubakov |J where it has been shown that in the large Nc limit, condensates in QCD are 
inhomogeneous and anisotropic at high fermion density. 

The Schwinger model is a favorite laboratory to test phenomena which are expected to 
happen in QCD^. In fact, an oscillatory inhomogeneous behavior in < ipip > was discussed 
in the Schwinger model |5| using operator bosonization. We think that the path-integral 
approach employed in the present paper is more appropriate to make apparent the crucial 
role that topological sectors play in the behavior of condensates (they are actually responsible 
for the non- vanishing of < x/ji/j >). In fact, our analysis implies that the phenomenon is not 
just a byproduct of 2-dimensional peculiarities. 

It is striking that the oscillatory behavior that we have found, exactly coincides (appart 
from the anisotropy that of course cannot be tested in one spacial dimension) with that 
described in || for QCD4. The structure of the N point correlation functions, given by 
eq . (f45[) , shows a non trivial dependence on spatial coordinates. This makes apparent that 
the ground state has, at finite density, an involved structure which is a superposition of 
standing waves with respect to the order parameter. 

Several interesting issues are open for further investigation using our approach. One can 
in particular study in a very simple way the behavior of condensates at finite temperature. 
The chiral anomaly is independent of temperature and plays a central role in the behavior 
of condensates through its connection with the index theorem. Therefore, one should expect 
(as discussed in ref. M for < tptp >) that a formula like ( fioD is valid also for T > 0. Of 
course, v.e.v.'s at — in the r.h.s. of this equation, should be replaced by those computed 
at finite temperature and hence the issue of zero- modes in a toroidal manifold should be 



carefully examined (see e.g. [p"6f ). 

Another extension which can be undertaken is related to the study of massless QCD 2 . 
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Indeed, the decoupling change of variables at the root of our approach, can be easily ex- 
tended for non-abelian gauge groups and has lead to deep insights in the properties of the 
model |J- 0. Finally, it should be worthwhile to consider massive fermions and compute 
fermion correlation functions, via a perturbation expansion in the fermion mass following 



the approach of [18]. We hope to report on these problems in a future work. 
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